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ABSTRACT 



Certain ORION missions may require three axis stabilization. Since ORION'S 
physical size severely limits its onboard fuel storage capability, passive stabilization 
techniques warrant investigation. This paper shows the development of linearized 
equations of motion and regions of stability with respect to gravity gradient stabiliza- 
tion. Gravity gradient stabilization by itself provides little yaw restoring torque: there- 
fore, additional torque generating devices are necessary to augment the gravity gradient 
effect. Control moment gyros, reaction wheels, and magnetic torquers will be investi- 
gated as to their suitability for ORION. 
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I. INTRODUCTION 



A. BACKGROUND ON ORION 

The philosophy behind the design and development of the Orion satellite is being 
driven in part by economic considerations. Small, affordable, general purpose satellites 
must be developed without sacrificing reliability. The preliminary groundwork on the 
feasibility of the ORION concept and design is detailed in [Ref. 1]. 

Orion is currently being designed to be launched from a Get Away Special (GAS) 
cannister aboard the Space Shuttle. The GAS cannister physically limits the size of 
ORION to approximately 35" in height and 19" in diameter. Figure 1 shows the Get- 
Away Special Cannister attached to the side of the shuttle bay. 




Figure 1. NASA Get-Anay-Special (GAS) Cannister [Ref. 1, p. 30] 



The Get-Away-Special (GAS) Cannister is designed to launch its satellite 
autonomously on command. A cross sectional view of the (GAS) Cannister is shown in 
Figure 2 on page 2. 

ORION should be able to meet the requirements of standard or typical payloads. 
Some of these requirements are listed by Reference 2 and given in Figure 3 on page 3. 
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Figure 2. Cross section of GAS Cannister [Ref. 1, p. 42J 



1. Background on attitude control 

The success of any satellite is dependent on its attitude control system. To be 
mission capable, the vehicle must first acquire a desired or specified orientation in space, 
and then maintain this orientaion within mission limits given such external factors as 
magnetic fields, aerodynamic drag, solar pressure, gravitational effects, and other dis- 
turbances. Deviations from desired orientations are detected through the use of sensors. 

Considerations for an attitude control system include: 
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(1) pointing accuracy 

(2) fuel consumption 

(3) thermal 



(4) power constraints 

(5) satellite lifetime 

(6) response time 



Mass 


3 2 lbm 


Volume 


2.36 ft 3 


Power 


34 watts 


DataRa te 


5000 bits/sec 


Orbit 


200-000 nm circular 


Inclination 


0 * -3 0 0 or 60' -120 ° 


Instruments 


Particle counter or Lens 



Figure 3. Summary of Ijpical Na\y/STP payload requirements 



The weight given to each of these considerations ultimately lies in the specific mission 
requirements designated for each satellite. The direction of this thesis will be to inves- 
tigate the feasibility of achieving sufficient performance utilizing a gravity gradient con- 
trol system. 1 he standard of +/- 1 degree for each axis will be used. This criteria is 
sufficient to enable ORION to satisfy most SIP (space test program) mission 
requirements.! Ref. 1, p. 93] 

do achieve stabilization, many methods of generating controlling torques have 
been developed. These methods are divided into either active or passive categories. Each 
of these categories has their relative advantages and disadvantages (Figure 4 on page 4 
and Figure 5 on page 4|. 

In this study, the primary focus is on a gravity gradient stabilization augmented 
by other torque generating devices to achieve the desired +/- 1 degree goal for each axis. 
A discussion of the equations of motion and gravity gradient stabilization theory is 
presented in Chapter 2. Disturbance torques and other torque generating devices are 
discussed in Chapter 3. Chapter 4 presents the results of stabilization schemes derived 
in Chapters 2 and 3. Observations and conclusions are discussed in Chapter 5. 
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Figure 4. Passive Stabilization Methods [Ref. 4, p. 19] 
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Figure 5. Active Stabilization Methods [Ref. 4, p. 19] 
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II. GRAVITY GRADIENT STABILIZATION THEORY 



A. COORDINATE DEFINITION 

Before beginning any discussion on attitude control a means of defining the coordi- 
nate system is necessary. The standardized coordinate axes for gravity gradient systems 
are as follows: the x axis points in the direction of satellite motion, the y axis is normal 
to the orbit plane, and the z axis points towards earth's gcocenter. These axes corre- 
spond to the roll, pitch and yaw axes, respectively. Figure 6 from Reference 3 defines 
roll (<}>), pitch ( 0 ), and yaw (>//) attitude errors as well as showing the standardized co- 
ordinate system. 
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13. COORDINATE TRANSFORMATION 



The relationship between two coordinate systems is essential for attitude control. 
There are methods available to describe this relationship of which directional cosines, 
quaternions, and Euler angles will be briefly discussed. Figure 7 on page 6 shows various 
methods of coordinate transformations along with their respective advantages and dis- 
advantages. 
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Figure 7. Coordinate Transformation Methods [Ref. 4, p. 412] 



1. Direction cosines method 

The directional cosine method involves a 3 x 3 rotation matrix A tJ where 



where 





A\\ A n 


V 




’a'.* 


X*y 


1 

N 

• 


^ij 


^21 ^22 


An 


= 


Y-x 




Y-z 




_A$ 1 ^32 


1 

I'D 

YD 




Z • x 


z-y 


Z • z 



X • x = cos 6 cos \p 

A' »y = cos ip sin 6 sin <p + sin ip cos <p 
X • z = sin \p sin 4> ~ cos \p sin 6 cos <p 

Y * x = — sin i p cos 6 

}' • j’ = cos (p cos \p — sin \p sin 6 sin <p 

Y • z = sin \p sin 6 cos <p + cos i p sin 0 
Z • x = sin 6 

Z • y = — cos 6 sin <p 
Z • z = cos 0 cos 6 



( 1 ) 
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2. Quaternions 

The quaternion method or as it is sometimes called, the Euler Synunetric Pa- 
rameter method, uses the Euler axis and rotation angle (a single rotation angle a about 

A 

a single axis £) method to describe the relationship between two coordinate systems 
[Figure 8 on page 7). 




The quaternion is a four component vector containing the same information as the Euler 
angle/axis transformation in which one coordinate axis is related to another and is of the 
form 



t hi , ft*, <743 



( 3 ) 



where 



r~ ■ a 

q x = E x sin — 
q 2 = E y sin — 

q 3 = E 2 sin — 

a 

q A = cos — 

The symbols i, j, and k satisfy the following conditions [Ref. 4, p. 758) 
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and 



.2 .2 ,2 , 

l = J =k = 1 

ij = -ji = k 
jk = -kj = i 
ki = -ik =j 



q] + q\ + q\ + q\ = 1 

When used for control purposes, (Ref 5] defines the first three elements of the 
quaternion {q u q 2 . q 2 ) as the respective roll, pitch, and yaw rotational errors: 

C troth huch' = 2 Ul' 72 < 73^1 

= sin sin , E 2 sin J (7) 

~[ 2 £’ JC , ty.Ey, o.E 2 ~\ 

The primary advantage of the quaternion method is that the computational time 
can be reduced by more that 40% over the equivalent directional cosine matrix solution 
in operations that require successive coordinate system rotations [Ref 5]. The transfor- 
mation from coordinate system A to coordinate system B can be described by the 
quaternion q = C<7d, qj- q^k. and the second transformation from system B to C can 
be described by q' = \_q\i. q'j< q\k, q\~\. The transformation from system A to system 
C may be described by the quaternion q" = \_q'\i, q"J, q " 3 , q" where 

q" = qq ' = (q { i + q J + q 3 k + q A ){q' x i + q'jj + q' 3 k + q\) ( 8 ) 

Multiplying out and substituting in the conditions stated above, the matrix form be- 
comes 



7”i 




q‘ 3 — ?2 


<■, fi 


7 1 


c i" 2 




7*3 4 A <k\ 


7' 2 


72 


<7"3 




7 ' 2 ~ < j\ q' a 


73 


73 


q'\_ 




-q'i ~q' 2 ~°i 3 


q ' a _ 


_ 74 _ 



Because of the quaternion property of interchangeability, a total of 16 multiplications 
and 12 additions are required to completely define the quaternion describing the trans- 
formation from coordinate system A to coordinate system C. To accomplish the same 
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transformation using directional cosines would require multiplying two 3x3 matrices, 
involving 27 multiplications and IS additions. The computational time savings for suc- 
cessive rotations transformations becomes apparent. The major drawback for the use 
of quaternions is that the numbers of the quaternion do not physically represent the 
transformation from one system to another. 

3. Euler angles method 

The Euler angles utilize three different rotations (angles), defining the orien- 
tation of a body with respect to an inertial reference frame. These angles are defined as 
<p (about the x axis), 6 (about the y axis), and ip (about the z axis). These axes, by es- 
tablished convention, may also be referred to as 1 (for x), 2 (for y), and 3 (for z) axes, 
respectively. 

The 321 Euler transformation described below starts with the inertial coordinate 
system XYZ. Figure 9 on page 10. a. shows the first rotation about the Z axis through 
an angle ip, which produces .v L v,Z. Figure 9 on page lO.b. shows the rotation of 6 about 
y, , producing .vy,z,. Figure 9 on page 10.c., showing the rotation of (p about x, produces 
xyz. The rotation described in Figure 9 on page lO.a. about the z axis produces 
Figure 10 on page 10. The transformation from XYZ to xy, Z is 

.Vj = X cos ip + F sin ip 
y, = — X sin \J/ + Y cos ip 

Z = Z 



or written in matrix form 



V 




cos ip 


sin ip 


o' 


~X~ 


y\ 


= 


— sin i p 


cos \p 


0 


Y 


z_ 




0 


0 


1 


_Z_ 
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Figure 9. Euler Angles 




Writing the transformations of Figures 9.b. and 9.c. in matrix form yields 



X 




cos 6 


0 


— sin 9 


~ x \~ 




y i 


= 


0 


1 


0 


y \ 


(12) 


_*i_ 




sin 6 


0 


cos 0 


_z_ 
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and 



„r 




"i 


0 


0 


X 


y 


= 


0 


COS 0 


sin 0 


y\ 


z _ 




0 


— sin 0 


COS 0 


j\_ 



( 13 ) 



Combining eqs( 11-13) into one equation yields the complete 321 coordinate transfor- 
mation from the inertial reference frame to the body reference frame 



~ X 




'i 


0 


0 


cos 6 


0 


— sin 6 


COS 0 


sin 0 


o" 




y 


= 


0 


COS 0 


sin 0 


0 
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— sin 0 


cos 0 
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J 


_2 _ 
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— sin 0 


COS 0 


sin d 


0 


cos 6 


0 


0 


1 


K_ 



cos 6 cos 0 

- cos 0 sin 0 + sin 0 sin 6 cos 0 
sin 0 sin 0 + cos 0 sin 6 cos 0 



cos 0 sin 0 

cos 0 cos 0 + sin 0 sin 0 sin 0 
— sin 0 cos 0 + cos 0 sin 0 sin 0 



(14) 
— sin 0 
sin 0 cos 0 
cos 0 cos 0 



A 7 

y 

z 



This is just one of the 12 different Euler coordinate transformations schemes available. 
Figure 11 on page 12 lists all possible Euler angle combinations and their resulting 
transformation matrices. It should be noted that eq(14) does not agree with the 321 
Euler transformation matrix of Figure 11 on page 12. The 321 matrix of Figure 11 on 
page 12 contains a error that interchanges 0 with 0 and 0 with 0. 
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CISIS# -SIC# 


-SI 


CIC# 


CIS# 




Slcl 


SISIC# -CIS# 


SISIS# * CIC# 


CIC# 


s# 


-CIS# 




-CISIC# » sisj 


CICI 


CISIS# * sic# 


s# sic# • c#s# 


-SIC# 


-S#S«S# ♦ CIC# 


CICI • $#S#S# 


SIC# 


-CIS# • SISIC# 


-SIC# • CISIS# 


cec# 


SIS# * CISIC# 


CIS# 


-SI 


CIC# 




CIC# -SISIS# 


CIS# • SISIC# 


-SIC# " 




-CIS# * 


CIC# 


SI 




SIC# » Cl MS# 


SIS# -CISIC# 


CICI 




CIC# 


CIS# 


-SI 




-CIS# » SISAC# 


CIC# • s#s#s# 


SIC# 




SIS# • CISIC# 


-SIC# • CISIS# 


CIC# 





MAtHIX A 



c» 

sis# 

CIS# 

c# 

-cis# 



SAS# 

CIC# -SICIS# 
-SIC# -CIC05# 

SIC# 
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•SIC# 

cis# * sicic# 
•SIS# • CICIC# 
SIS 4 

CICIS# * SIC# 



SIS# 


-SICIC# -CIS# 


-sicis# *c#c# . 


CIC# -SICI5# 


SIS# 


-cis# -sicic# 


SOS# 


c# 


soc# 




Sir# « CICIS# 


-CIS# 


-SIS#* CICIC# . 


CICIC# -SIS# 


CIS# 


-CICIS# - 


SIC# 


-SIC# 


c# 


SIS# 




SICIC# « CIS# 


SIS# 


-SICIS# * 


CIC# . 


CIC# -SICIS# 


CIS! ♦ SICIC# 


SIS# ' 




-SIC# -CICIS# 


-SIS# • CICIC# 


CIS# 




SIS# 


-SIC# 


c# 




CICIC# -sis# 


CICIS# ♦ SIC# 


-CIS# 




-side# -cis# 


-S#CIS#*C|C# 


IIS# 




sic# 


SIS# 


Cl 





Tlguie tl. '1 able of Uuler angle combinations IRef. 4, p. 764] 



12 



C. RIGID BODY EQUATIONS 

The general motion of a particle and its corresponding momentum have been written 
about extensively in various publications. A summary of the relevent derivations and 
equations of motion from [Ref 6] are presented here. 

1. Derivative of a Vector in a Rotating Coordinate System 

The time derivative of a vector in a rotating coordinate system consists of the 
rate of change of the vector relative to the moving axis and the rate of change due to the 
rotation of the axis. Figure 12 shows a vector in a rotating system where XYZ represent 
an inertial system and xyz are the set of axes rotating with angular velocity a> relative 
to XYZ. If unit vectors along the X, Y, and Z axes are i,j, and k, respectively, the vector 
r can be written in the form 



r = xi + yj + zk (15) 




The time derivative of r is 



r — xi + xi + yj + yj + zk + zk 



(16) 



From [Ref. 7j, it has been shown that the time derivatives of the unit vectors are 



1 = (O x I 



(17) 
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(18) 



j = w xj 

k = u) x k (19) 

Defining r„, as the time rate of change of r relative to xyz, eq(16) becomes 

r = r rel 4- cl> x r (20) 

2. General Particle Motion 

Translational motion of a particle with references to two coordinate systems is 
shown in Figure 13. 




Figure 13. Translational motion of a particle [Ref. 6, p. 107] 

The position vector of a particle relative to coordinate system XYZ is 

R = R c + r (21) 

where R c is the vector from the axis of the XYZ coordinate frame to the xyz axis frame 
and r is the differential (position) vector in the xyz frame. The velocity of the particle 
is the derivative of eq(21) with respect to time 

V=R = R c + r (22) 
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The velocity at the origin R c = V c combined with eq(20) yields 



V=V c +V rel + coxr (23) 

3. Momentum of a Rigid Body 

A particle's momentum is the product of its velocity and mass. The linear mo- 
mentum of the particle of Figure 13 on page 14 is 

P = I V dm (24) 

where V is the particle's velocity and m is its mass. Recalling that V re , is the time rate 
of change of r relative to the xyz coordinate system, and assuming that the xyz axis is 
centered in the rigid body, V re/ = 0. Substituting eq(23) into eq(24) with F, e/ =0, 



P = I (K c + wxr) dm = mV c + oj x I r dm (25) 

^ m ^ m 

By choosing the center of mass of the rigid body as the center of the rotating coordinate 
system, / r dm = 0 by definition. This simplifies eq(25) to 

P — m V c (26) 

where P is the linear momentum of the rigid body, m is its mass, and V c is the velocity 
of its center of mass. 

The angular momentum of a particle dm about point C is defined as the moment 
about C of the linear momentum given by 

h c = r x V dm (27) 



Substituting for V yields 



H c = I r x V dm = I r x{V c + (o x r) dm 



(28) 



Rearranging terms 
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(29) 




r dm x 



V c + 



[ 



r x (co x r) dm 



By choosing the center of mass of the rigid body as the origin of the rotating coordinate 
system, eq(29) reduces to 



H c = r x (co x r) dm (30) 

^ m 

Reducing the angular velocity vector co into its components 

co = co x i + coy + co.A' (3 1 ) 

Combining equations(15), (30), and (31) yields 

He = Uxx°>x - IxyOy ~ 

+ £ l X y(0 X *b iyyCOy -CO^^/ (32) 

+ [ l X M X iy-COy +I,.(0 Z ']k 

where 



I XX 



(/ + z 2 ) dm 



(33) 



lyy 



= j (- V2+Z2 ) 



dm 



(34) 



*22 



-JV+/) 



dm 



(35) 



are the moments of inertia about the x,v, and z axes, respectively, and 



‘xy 



= J xy dm 



(36) 
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(37) 



-I* 



I X2 = | xz dm 



=L 



lyz = I vz dm 



(38) 



are the products of inertia. Equation(32) written in matrix form becomes 



1 




H y 


= 


1 

1 





I xx ^xy 

—*xy lyy 

~^x 2 ~ ly: 



~h 

-l 



yz 



L 



a). 



a>„ 



L w ;. 



(39) 



4. Equations of motion/Euler's Moment Equations 

The derivation of the following Euler's moment equations are summarized from 
[Ref 6). Recalling that angular momentum h c relative to the point C is 



= r x mv 



(40) 



its derivative with respect to time is 



h c = r x mv + r x mv (41) 

The moment of the force about a point C by definition is 

M c = r x F (42) 

Recall Newton's second law describing force: 

F = ma = mv (43) 

Substituting eqs(22) and (43) into eq(41) yields 

h c = r x m(R c + r) + r x F (44) 

Substituting eq(42) into eq(44) and rearranging terms yields the moment equation 
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M C = H C + 



(R c x /•) dm 



(45) 



where H c is the angular momentum of the rigid body about point C. From eq(20) the 
latter part of eq(45) can be written as 



R c x r dm = R c x 




(46) 



Noting that point C is the center of mass f r dm = 0 and r re , = 0 for a rigid body, eq(46) 

* m 

reduces to 



U c =il c (47) 

The moment of external forces about the center of mass of a rigid body is equal to the 
time rate of change of the angular momentum of the body about the center of mass [Ref. 
7]. In terms of a rotating axis, the moment equation of eq(47) becomes 



M c =H crel +(oxH c (4S) 

Substituting eq(39) in for H, the components of the moment equation can be written as 

M x =H x + co y H 2 -co 2 H y (49) 

M, = il y + co 2 II x - CJ x H 2 (50) 

M 2 = //. + 0J x II y — co y II x (51) 



where are the moment equations. In this study, an assumption is made that the 
principle axes of the moment of inertia are the rotating axes x,y,z; therefore, the pro- 
ducts of inertia (I xy , I y2 ) are equal to zero. The angular momentum components re- 
duce to 

U.x = (52) 
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( 53 ) 



Hy = iyyXOy 



U 2 = I. z 0J. 



(54) 



Substituting eqs(52-54) into the moment equations yields Euler's moment equations 

M X = ! XX^X + OJyOJ Z (I ZZ - Iyy ) (55) 

My — IyyiOy H" 10^(0 -( if.) (56) 

47- == / 77 CU 7 + (l) X (Oy{Iyy 1 X X' ( 5 ^) 

Referring to the coordinate system shown in Figure 6 on page 5, the angular 
velocities in terms of orbital rate co 0 and attitude error angles pitch(0), roll(0), and yaw 
(0) from [Ref. 6, p. 130] are 



(Or 



(O v 



L^zJ 



1 0 — sin 9 0 

0 cos 0 cos 6 sin 0 6 

0 — sin 0 cos 6 cos 0 0 

cos 6 sin 0 

— a> 0 cos 0 cos 0 + sin 0 sin 6 sin 0 
— sin 0 cos 0 + cos 0 sin d sin 0 



When using small angle approximations, eq(5S) reduces to 



(58) 



°>x 




1 

0 

3 

1 

-e* 

i 


iOy 


= 


0 -(O 0 






0 + O>0 <t> 



(59) 



Note that the orbital rate is defined as 




(60) 



where n is the gravitational constant and R is the radius of the earth plus the satellite's 
altitude. For the purpose of station keeping, the small angle approximations are valid. 
Substituting eq(59) into the rigid body moment eqs(55-57) yields the linearized rigid 
body dynamic equations 
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(61) 



^ xRB I.xxfi d~ ^0 (/>>■ 7-..)0 ^oi^xx "1" lyy)*^ 



MyRB ~ ! yy e 



(62) 



•• 2 • 

^iRB = ^11^ "1" < ^*o(A’>' ^xx)^ ^o(^XX "1" ^ II Iyy)4* 



(63) 



D. GRAVITY GRADIENT EFFECT 

Attitude control using gravity gradient stabilization works on the principle that 
given an object with an asymmetrically distributed mass, the object will tend to align it- 
self with its minimum moment of inertia axis along the local vertical towords the earth. 
The force due to gravity on a mass element is 



where G = universal gravity constant, M e = the mass of the earth, g e = gravity acceler- 
ation at the earth's surface, and R e = the radius of the earth. Reference 8 shows the 
torque generated about the mass center of the object is 




(64) 



where 



R = Rq + r 



(65) 



and 



ii e = GM e = g e R ; 



(66) 



dM g = r x dF = — 



Hedm 

— rxR 

R 2 




(67) 



Cubing eq(65) yields 




3. rR, 3 SRa r * 




(68) 
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Noting that r<^R 0 , the higher order terms drop out. Inversion of eq(68) yields 



R Rn 



♦ R 0 
Ro 



( 69 ) 



Substituting eq(69) into eq(67) and integrating to find the moment yields 



f f 

M g = — n e \ r dm x R 0 H — 

Rn 



3n e f (r x R 0 )(r • R 0 ) 



o m 



Ro 



dm 



(70) 



Given that the center of mass is the origin of the body frame, f r dm = 0, eq(70) reduces 
to 



Mg = f(rx R 0 )(r . R 0 ) dm 
Rn J m 



(71) 



In the reference system defined by Figure 14, the K (Z) vector is pointing towards 
the earth; therefore, the vector R 0 in the reference coordinate frame is 



R 0r =RoK (72) 




Figure 14. Gravity gradient torque [Ref. 6, p. 113J 
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Restating the Euler 321 rotation matrix described earlier 



/ 




"l 


0 


0 


cos 0 


0 


— sin 9 


COS \p 


sin \p 


o' 




7 


= 


0 


cos 4> 


sin 4> 


0 


1 


0 


— sin \j/ 


cos 4/ 


0 


J 


_A-_ 




0 


— sin 4> 


cos 4> 


sin 0 


0 


cos 9 


0 


0 


1 


_A_ 



COS 0 COS 


cos 9 sin 4 / 


(73) 
— sin 0 


"/' 


— cos 4> sin 4> + sin 4> sin 9 cos 4> 


cos 4> cos 4 1 + sin </> sin 9 sin 4* 


sin 4> cos 0 


J 


sin 4> sin 4> + cos 4> sin 9 cos i p 


— sin <p cos 4* + cos </> sin 0 sin 4* 


cos 4> cos 0 


A 



Transforming R 0r into the spacecraft's fixed coordinates i. j, and k: 





COS 9 cos 4* 


cos 9 sin *p 


— sin 9 


'o' 


— 


— cos 4> sin 4> + sin <p sin 9 cos \p 


cos 4> cos 4* + sin 4> sin 9 sin ip 


sin 4> cos 9 


0 




sin 4> sin 4* + cos 4> sin 9 cos ip 


— sin 4> cos 4> + cos 4> sin 9 sin i p 


cos 4> cos 9 


_ A”_ 



which reduces to 

R 0r = R 0 ( — sin 9i + sin 4> cos 9j + cos 4> cos 9k) (75) 



Substituting eq(75) into eq(71) yields the gravity gradient moment equations 



M„ = 



W e 

Rn 



(A, — I yv ) sin 4> cos 4> cos ‘0 

— (I xx — I, z ) sin 9 cos 9 cos 4 > 

— (Iyy — I ., w ) sin 9 cos 9 sin 4> 



(76) 



Replacing 



2 



(77) 



where cu 0 is the orbital rate, and making the small angle approximations, the gravity 
gradient moment equations become 



M g = 3coq 



( kyy ) (P 

(/»-/«)* 



0 



(78) 
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Equation(78) clearly shows that by increasing the moments of inertia for I x and l y 
with respect to 1„ the gravity gradient moments for roll and pitch will increase. By ex- 
tending a point mass along the z axis, a sizeable increase in the moments of inertia along 
the x and y axes is realized. 

The magnitudes of the gravity gradient torques generated at different altitudes is 



given by Figure 15. 




E. DESIGN FOR STABILITY 

1. Stability and Frequency Analysis 

Combining the rigid body dynamic equations (61-63) with the gravity gradient 
equation (78) yields the undamped linear dynamic equations for roll, pitch, and yaw: 

Roll I x ip + 4c o 2 0 (I y - 1 2 )4> - oj 0 (I x +I 2 - I y )j, = 0 (79) 



Pitch lyd + 3cuq {I x — I 2 )6 = 0 



(80) 



23 



( 81 ) 



Yaw I 2 ij/ + a)\(l y - l x )ip + oj 0 {I x + I 2 - l y )4> = 0 
Using Reference 9 notations for simplicity, the following variables are defined: 

I y ~ h 

A=—r— (82) 

J X 

I x -L 

B = (83) 

‘y 

Iy ~ I x 

C = JL 1 — (84) 

Substituting eqs(S2-S4) into eqs(79-Sl) yields the following second order differential 
equations 



0 + 4/Jo»o0 — (1 — A)a) 0 i// = 0 



(85) 



6 + 3Bco 0 d — 0 



(86) 



0 + Cto o 0 + (1 — 00 = 0 



(87) 



Taking the Laplace transform of eqs( 85-87) allows for frequency analysis. The results 
are put into matrix form: 



S~ + 4.4 coq 0 


-(1 ~A)co 0 S 


V 




"o“ 


0 S 2 + 35a»o 


0 


e 


= 


0 


(l-C)a>oS 0 


1 

CM O 

a 

+ 

rs 


A 




_o_ 



( 88 ) 



The determinant of the matrix is the system's characteristic equation. Using the middle 
row’ to simplify calculations, the characteristic equation is found to be 

(S 4 + ( 1 + 3 A + CA)coqS 2 + 4 A Ccoq) (S 2 + 3 Bcoq) = 0 (89) 

"Solution of the first term on the left yields the natural frequencies of roll and yaw, 
repectivelv; the second term on the left defines the natural frequency of the pitch axis." 
[Ref. 9, p. 9] Solving for the first term on the left (roll/yaw stability) by the quadratic 
formula and dividing through by a>^ yields 
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( 90 ) 



SL S yaw ~ (1 + 3 A + AC) ± N /(l + 3A + ~ACj - 1 6AC 

2 2 2 

COq COq 

For roll/vaw stability the roots of eq(90) must be real and negative; therefore the fol- 
lowing stability conditions exist [Ref. 4 p. 611] 

1 +7>A + AC>4j~AC (91) 

AC> 0 

Solving for the second term of Eq(89) and dividing through by <yj] yields 

S pitch „ n 3(C— A) 

— = 

The condition for stability for the pitch axis is 

A>C (94) 

Converting A, B, and C back into their respective moment of inertia relationships de- 
fined by eqs(82-84), combined with the conditions for stability, yields two possible ori- 
entations: I y > I x > I 2 and I x > I 2 > I y . Figure 16 on page 26 shows a summary of the 
stability regions based on moment of inertia ratios. [Ref. 10 p. 204] 



(92) 



(93) 
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Figure 16. Gravity Gradient Stability Regions 
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III. ENVIRONMENTAL DISTURBANCES AND RESTORING TORQUES 



A. ENVIRONMENTAL DISTURBANCES/TORQUES 

The are four major environmental disturbances that affect an orbiting spacecraft: 
aerodynamic, magnetic, gravity gradient, and solar radiation torque. Figure 17 lists 
these disturbances and their relative regions of dominance. 



Source 


DEPENDENCE ON DISTANCE 
FROM EARTH 


REGION OF SPACE WHERE 
OOMINANT • 


AEROOYNAMIC 


• 


AITITUOES BELOW - 500 km 


MAGNETIC 


I/, 3 


1 - 500 km TO - 35.000 km; 1 


gravity Gn aoient 


I/. 3 


1 II E..OUT TO ABOUT SYNCHRONOUS ALT I T UDE 1 1 


SOLAn RADIATION 


INDEPENDENT 


interplanetary spaqe above synchronous 

ALTITUDE 


MICnOMETEoniTES 


LARGELY INDEPENDENT; HIGH CONCENTRATION 
IN SOME REGIONS OP THE SOLAR SYSTEM 


NORMALLY NEGLIGIBLE; MAY BE IMPORTANT IN 
SOME SMALL REGIONS lINTERlon OF SATUnN S 
RINGS! 


ALTITUDES LISTED ARE ONLY REPRESENTATIVE; THE SPECIFIC ALTITUDES AT WHICH VARIOUS TOnOUES DOMINATE ARE HIGHLY 

spacecraft dependent. 



Figure 17. Environmental Disturbance Torques [Ref. 4, p. 17) 



1. Aerodynamic disturbance 
a. Drag 

"For satellites near the earth the principle non-gravitational force is aero- 
dynamic drag. Aerodynamic drag is a retarding force due to atmospheric friction and is 
in the direction opposite the space vehicles velocity vector." [Ref. 4, p. 63] 

The drag force is a function of vehicle velocity, air density, and surface area 
of the satellite. The drag force is given by the following equation: 



pv'Cja 

2 



(95) 



where d is the drag force, p is the atmospheric density, v is the velocity of the satellite 
along its orbital path, c d is the drag coefficient, and a is the area of the satellite over 
which the drag acts. 

b. Aerodynamic torque 

The collisions of air molecules of the upper atmosphere with the satellite 
surface produces a torque about its center of mass. Aerodynamic disturbances are a 
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function of the satellite's altitude, velocity, and symmetry. A satellite with a relatively 
low earth orbit may be significantly affected by aerodynamic torques. 



center of pressure (cp). The offset between the center of mass (cm) and the center of 
pressure (cp) is denoted by L x , /,„ and L„ respectively. The simplified aerodynamic dis- 
turbance torques from Ref 8 may be computed by the following expression: 



where T„ is the torque due to aerodynamic pressure (ft-lbs), P„ is the aerodynamic pres- 
sure (lb sift 1 ), l„ is the distance between the center of mass and the center of pressure (ft), 
A, is the exposed surface area (ft 2 ) and a is the angle of attack (radians). For practical 
purposes, the aerodynamic torque in the x direction is zero (sin a = 0 when asiO). The 
conversion factor from ft-lbs to N-m is approximately 1.356. Figure 19 on page 29 
shows aerodynamic pressure as a function of altitude. 



Shown in Figure 18 is a cylindrical satellite with its center of mass (cm) and 



T a = PJa A a sin a 



(96) 




R 




Z 



Figure 18. Aerodynamic effects for cylindrical satellite [Ref. 3, p. 58] 
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Figure 19. Aerodynamic pressure as a function of altitude [Ref. 8, p. 456] 



The magnitude of the pitch and yaw aerodynamic disturbances for a one 
inch cm-cp offset are shown in Figure 20 on page 30 for various altitudes. The magni- 
tude of these torques decrease significantly with increases in altitude. "Proper control of 
the satellites inertias to enhance gravity gradient torques can make aerodynamic dis- 
turbances essentially insignificant for orbital altitudes above 300 nautical miles (480km). 
On the other hand, in the 100 to 200 nautical mile (160-320km) altitudes aerodynamic 
disturbances are 3 to 4 orders of magnitude greater than the gravity gradient torque from 
a 1 slug-fr J inertia difference and 1 degree attitude error. In the lower regions, very large 
inertia differences are required for accuracy. [This leads to design considerations] where 
proper configuration can make these aerodynamic torques work as restoring torques 
instead of disturbance torques" [Ref. 3, p. 57]. 
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Figure 20. Aerodynamic torque for a 1 inch offset as a function of altitude. 

2. Rigid Body Solar Pressure Torque 

"The major factors determining the radiation torque on a spacecraft are the in- 
tensity and spectral distribution of the incident radiation, the geometry of the surface 
and its optical properties, and the orientation of the sun vector relative to the spacecraft" 
[Ref. 4, p. 570]. 

Solar pressure torques are the disturbances produced by solar radiation pres- 
sure and is a function of the offset of the center of pressure (cp) and the center of mass 
(cm). Referring to Figure 21 on page 31, the center of mass-center of pressure offsets 
are designated as L„ L y , and L 2 , respectively. The torques generated due to solar pres- 
sure are given by [Ref. 8, p. 454| as: 

T s = P S L S A S cos e (97) 

T s = 2P S L S A S cos 2 e (98) 

where T, = torque due to solar radiation (ft-lbs), P, = radiation pressure (lb/f/ J ) and is 
s; constant at 9.65 x 10* 8 lb ,'J't 1 for an earth orbiting vehicle, L, = center of mass center 
of pressure offset (ft), A, = surface area of satellite normal to sun(/P), and t = angle 
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of incidence (degrees). Note, equation (97) is for an absorbent satellite body and 
equation (9S) is for a reflective body. 




Figure 21. Rigid body solar pressure torques [Ref. 3, p. 61) 



B. RESTORING TORQUES 

As seen from chapter 2, the restoring torques of the undamped rigid body combined 
with the gravity gradient effects are considerably smaller for the yaw axis than for either 
the pitch or roll axes. This is due to the necessity of large moments of inertia along the 
pitch and roll axes, as compared to the yaw axis, for gravity gradient stabilization to 
occur. The lack of adequate yaw restoring torque requires that additional yaw restoring 
torque be provided for accurate 3-axis stabilization. A discussion of yaw restoring torque 
methods follows. 

1. Thrusters 

Thrusters are an effective torque generating device; however, for the purpose 
of this thesis, thrusters will not be considered for attitude control. Thrusters are neces- 
sary for orbit maintenance and orbital transfer. Due to the physical size constraints im- 
posed by the GAS cannister, ORION is severely limited in its fuel storage capacity to 
approximately 71 lbs of hydrazine fuel [Ref 1]. An assumption that the launch platform 
will place ORION in its final orbit is optimistic. Figure 22 on page 32 shows the fuel 
required to go from one orbital altitude to another. Saving the onboard fuel for orbit 
maintenance and transfer and not for attitude control allows more fuel to be devoted 
towards supporting the mission either by allowing greater orbital flexibility and/or a 
longer lifetime. 
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Figure 22. Fuel required for orbit transfer based on a 2501b satellite. [Ref 2, p. 23J 
2. Momentum Exchange Devices 

The principle behind momentum exchange devices is that by changing the an- 
gular momentum of the control device, the angular momentum of the vehicle will be 
changed an equal and opposite amount. Recall that the rate of change of angular mo- 
mentum is equal to the torque that is generated. The two basic types of momentum ex- 
change devices are reaction wheels and control moment gyros (CMG). The basic 
difference between the two is the method by which the change in angular momentum is 
accomplished. The reaction wheel axis is fixed (generally along one of the axes of the 
vehicle) so that the change in angular momentum is accomplished by varying the speed 
of rotation of the flywheel. The CMG creates the change in angular momentum by tilt- 
ing the contain speed flywheel with respect to the vehicle. 
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a. Reaction wheels 



"An axis may be controlled by varying the speed of the reaction wheel in 
response to an attitude error [Ref. 6, p. 1 49). The momentum of a system with a reaction 
wheel is 



H = H W +H C 



(99) 



where H c is the angular momentum of the system and 



H w - 






wy 



hL 



is the angluar momentum of the reaction wheel. Aligning the reaction wheels along each 
axis, and substituting eq(99) into eq(4S) and adding the gravity gradient moments of 
eq(7S) yields 



I x 4> + A<o\{I y - I 2 )<f> - (o 0 {I x + I. - I y )\p + (6 - (o 0 )H wz - (i '(/ + to 0 (j))H uy + H wx = ^T x 

I y 6 + 3fc>o U x — I Z )Q + ( — <t> + co 0 t/')//*., + (</' + co 0 4>)H wx + H wy = ' S ^_T y (101) 

I 2^ + WoVy ~ 4)'A + + I z — I y )4> + (4> — co 0 i I/) + ( — <j> + co 0 ip)H wx + H W2 = /^T. 

"Because H„ x ,H wy ,H W2 , and co 0 are small, the coupling terms are small. If the coupling 
terms are neglected, the equations of motion about the roll, pitch, and yaw axes become 
independent, and hence, they can be controlled independently" [Ref. 6, pp. 149-150]. 
Dropping the coupled terms, the motion equations become 

I x 4> + 4a )\{I y — I 2 )4> — c o 0 {I x + I 2 — I y ,)<A 4- H wx = '^T x (102) 

I y e + 3 a>l(I x - I 2 )6 + H wy = Y,Ty (103) 

W + co 2 0 (Iy - I x )iP + co 0 (I x +I 2 - Iy)4> + H W2 = Y/z (104) 
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Allowing the momentum of the reaction wheel to be proportional to the attitude error, 
then by differentiation, the reaction wheel control torques are equal to 



J xrw = ~ K x (t x (}> + 4>) 



(105) 



Tyrw = — ky^Tyd + 6 ) 



(106) 



7 2TW 4* *A) 



(107) 



and the control scheme would be as shown in Figure 23. 




Figure 23. 3 axis reaction wheel control scheme [Ref. 6, p. 149] 



The torque genterated by the reaction wheel is a function of the wheel's angular velocity 
and its moment of inertia. This means that a small lightweight wheel at a high rate of 
speed generates the same amount of torque that a large wheel at a low rate of speed. 
Since ORION is designed for a relatively short lifespan, wear and tear on spinning parts 
is not an overriding concern. It should be noted that continuous disturbance torques 
that are not cyclical in nature will eventually saturate the reaction wheel, whereby, the 
reaction wheel rotational speed will reach a maximum. Thrusters or some other method 
of momentum removal must be used to reduce wheel speed. Disturbance torques that 
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are cyclical in nature arc easily offset by the reaction wheel as long as the disturbance 
torque is such that it does not saturate the reaction wheel speed. 

3. CMG(s) 

"A control moment gyro's (either single or double gimbaled) angular mo- 
mentum is due to the rotor which is spinning about the spin axis with a constant 
angular rate. Because the spin axis is gimbaled, a commanded gimbal rotation causes 
the direction of the angular momentum vector to change, thus creating a control 
torque parallel to the output axis. The magnitude of this torque depends on the 
speed of the rotor and the gimbal rotation rate." IRef. 4, p. 200] 



A control moment gyro scheme is shown in Figure 24. 




Because of their expense and weight, CMG's are only used on large spacecraft 
[Ref. 4]. For this reason, control moment gyros will not be addressed. 

4. Magnetic torquers 

Magnetic torquers may be used to dampen the kinetic energy of the libration 
motion as well as correcting attitude errors in any direction. The magnetic torquers 
consist of three wire coil sets arranged perpendicularly to each other. This configuration 
allows magnetic torque to be generated along any axis or in any direction as a result of 
adding two or three different orthogonal torques. As shown in Figure 25 on page 36, 
when a current is applied around the loop of a coil, it produces a magnetic dipole that 
is normal to the plane of the coil and a magnitude which is proportional to the coil's 
enclosed area and ampere-turns. 
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The magnetic moment M is defined as 

Af = MAn (10S) 

where N is the the number of turns around the coil, 1 is the current in amperes, A is the 
area enclosed by the loop, and n is the unit vector normal to the plane of the loop. Be- 
cause of the relationship between N. 1, and A, many different combinations are available 
to produce the desired magnitude of the magnetic moment. Figure 26 on page 37 shows 
three existing spacecraft with various combinations of number of turns and current 
drawn. The torque generated by the magnetic torquer is 

r mag = MxB (109) 

where M is the magnetic dipole and B is the earth's magnetic field. The earth's magnetic 
field for an equatorial orbit can be approximately modeled as a tilted dipole located at 
78.3° north latitude and 69° west longitude as shown in Figure 27 on page 37. 
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SPACECRAFT 


REMARKS 

. ----- 


SAS-3 


COIL CONSISTS OF 260 TURNS. MAXIMUM 
CURRENT IS 0.6 A AND MAXIMUM POWER 
CONSUMPTION IS 10 W 


OSO— 8 


COIL CONSISTS OF 360 TURNS. MAXIMUM 
CURRENT IS 0.075 A 


AE — 3 


TWO COILS; EACH HAS 500 TURNS. MAXI- 
MUM POWER CONSUMPTION IS 12 W 



Figure 26. Some existing magnetic moment schemes [Ref. 3, p. 205] 
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Figure 27. Earth's magnetic field modeled as a tilted dipole [Ref. 6, p. 148] 



The strength of the earth's magnetic field for a circular equatorial orbit in component 
form can be approximated by 



37 



( 110 ) 



B x — sin(l 1.7) sin(281 — lam s ) 

R 

Mo 

B v f- cos(l 1.7) (111) 

R 3 

Mo 

B z = — f- sin( 11.7) cos(2Sl — larn s ) (1 12) 

R 

where lam s is measured in degrees of right ascension from Grenich, M t is the strength 
at the earth's center which is approximately 8 x 10 ls Wb-m, and R is the distance in 
meters from the earth's center to the orbital altitude. With respect to the satellite, the 
earth's magnetic field is changing continuously as the vehicle moves along its orbital 
path. Putting into component form both \1 and B, eq(109) becomes 

T mag - X = M : B y - M y B : (113) 

T mag _ y = M x B z -M z B x (114) 

T mas -: = M y B x - M x B y (115) 

It should be noted that the earth's magnetic field is not uniform and varies in 
direction based on the satellite's position. For accurate attitude control, the vehicle must 
be able to exactly know both its position and the earth's magnetic field at that position. 
For low earth orbits, the magnetic field is well mapped and may be easily programmed 
into the attitude control computer. For the purpose of this thesis, magnetic field ap- 
proximations will suffice to demonstrate magnetic torquing as a viable attitude control 
capability. 
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IV. CHAPTER 4 



A. GRAVITY GRADIENT STABILIZATION RESULTS 
1. Attitude Control 

There are three general attitude control phases with respect to gravity gradient 
systems: gravity gradient capture, transition period, and steady state. Gravity gradient 
from a tumbling mode occurs when angular velocity slows to one rotation per orbit or 
less, allowing the minimum moment of inertia axis to align with the local vertical. The 
transition period is the time from gravity gradient capture until a steady state has been 
achieved. This period of oscillations will vary depending on the dampening scheme used 
to reduce the amplitude of oscillations. Predominant throughout the satellite's lifetime, 
the steady state phase is the period in which the satellite's errors are confined within 
operationally acceptable amplitudes. Maintenance of the desired steady state requires 
minimization of external disturbances and control of those which cannot be eliminated. 

The following group of simulations assumes that a steady state phase has been 
achieved and that the small angle approximations are in effect. For each of the simu- 
lations. the following initial conditions are in effect: </>, 6 , and i p position error = 5°, and 
</>, 6 , and p = 0. 



2. Various results based on different configurations 

a. Effects of extending booms 

For Figure 28 on page 40. a^lO kg point mass is used and the mass of the 
boom is assumed to be negligible. Extension of the boom increases the moment of in- 
ertia in the x and y directions. This increases the restoring torque along the roll and pitch 
axes although it does little to dampen the libration motion. 

b. Gravity gradient effect alone 

The following figures show the gravity gradient effect on a satellite at an 
altitude of 1000 km having moments of inertias as shown. As shown by Figure 29 on 
page 41 and Figure 30 on page 42, stabilization along the roll and pitch axes would take 
a prohibitively long time while the yaw direction would not stabilize. It is clear that 
gravity gradient stabilization alone is incapable of acceptable three axis stabilization. 
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MOMENT OF INERTIA (KC/M"2) 



MOMENT OF INERTIA VS DOOM LENGTH 




Figure 28. Effect of extending boom on moments of inertia 
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Figure 29. Gravity gradient effect only Ix= 15, Iy= 18, Iz= 3 



41 




c. Effects of Reaction ))'7ieel(s) 

One method of augmenting gravity gradient stabilization is through the use 
of reaction wheels. Teldrix reaction wheel classes and corresponding technical data are 
presented in Figure 31 on page 43. 
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Wheel diameter 


cm 


20 


26 


35 


50 


60*) 


Angulai momentum range 


Nrns 


1 .0...6.5 


5.0... 20 


14. ..BO 


50.. .300 


200... 1000 


Max. reaction torque 


Nm 


0.2 


0.2 


0.2 


0.3 


0.3-0. 6 


Speed***) 


min r 


6000 


6000 


6000 


6000 


6000 


Loss torque at max. speed**) 
Power consumption: 


Ntn 


< 0 012 


< 0.013 


<0.015 


< 0.022 


< 0.07 


— steady stale (depending 
on speed) 


W 


2. ..7 


2...G 


2... 10 


3... 15 


10. .50 


— max. power rating 


W 


< 60 


< 00 


< 100 


< 150 


< 500 


Dimensions: 
- diameter A 


mm 


203 


260 


350 


500 


600 


- height B 


mm 


75 


85 


120 


150 


180 


Weight 


k 9 


2.7. ..3.4 


3.5. ..6.0 


5.0...B.0 


7.5... 12 


20. .37 


Environmental conditions: 














- operating temperature 

- vibration (sinusoidal) 

- vibration (landom) 
-linear acceleration 




suitable for satellites compatible with latrncheis 
such as ARIANE or Space Shuttle 



Figure 31. Teldrlx Reaction Wheel Teclinical Data 



Of the five diameter classes, the 20 cm diameter wheel is the most appropriate of the 
classes listed for ORION. Reaction wheels also come in diameters smaller than 20 cm 
diameter. Bendix Aerospace SP reaction wheel of Figure 32 on page 44 has a diameter 
of 6 inches. 

As illustrated in Figure 31 and Figure 32 on page 44, power consumption 
values dilfcr with different reaction wheels. For the 20 cm Teldrix reaction wheel, steady 
state requires 2-7 watts with a maximum less than 60 watts. The Bendix SP reaction 
wheel requires 2.5 watts at steady state with a maximum power consumption of 6.5 
watts. Due to weight and power consumption values, the Bendix reaction wheel appears 
to be the superior choice of the two. 

The following figures were generated by setting the reaction wheel torques 
equal to eqs( 107-109). The initial conditions are 0 , 6, and \J/ position errors = 5°, and 
4> = 6 = \jj = 0. For simulation purposes, the maximum reaction torque value of .2 N-m 
from Figure 31 will be in effect. 

Figures 33-35 have moment of inertias /, = 15, I y = 18, and /, = 3, with 
K, = K y = 2.0, K. \ = 1.0, and x x — x y = x,= 1.0. Figure 33 on page 45 shows <f> (position 
error), </> (velocity), and the restoring torque ( T rwi ) generated by the reaction wheel with 
respect to time. Figures 34 and 35 show 6 , 8 , T nt , and \j/ , T rw ^, respectively. With 
the relatively small moments of inertias, stabilization of all axes within the ± 1° criteria 
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occurs within 40 seconds. Increasing the moments of inertias to / 2 = 95, 7^ = 98 , and 
7 2 = 3 while holding K„ K y at 2.0, K, at 1.0, and t 2 , x y , and r 2 constant at 1.0, yields 
Figures 36-3S. As expected, the 0 and 0 axes take longer to stabilized (approximately 
160 seconds) than the above case with the smaller moments of inertia. Th£ 0 axis of 
Figure 38 stabilizes much faster than either the 0 or 0 axes. This is due to the relatively 
small moment of inertia (7 2 ) compared to the restoring torque generated by the reaction 
wheel. 



rROGRAM/VNIT NAME : 


SP Reaction Whs • 1 


UM1T TYPB MO. (ARf) | 


3090003 


MOMENTUM ( FT -LR-SBC ) i 


Q.S 1 1250 RPM 


V B l CUT (1.DS.) j 


*.9 


SUR : 


6.0" D1A. MAX 1 3.9" HIGH (excluding counting 
feet projection) 


STALL torque ( OZ - 1 M ) 1 


2.2 


r EAX POWER (WATTS) : 


6.5 


MAX STEED ( RPM ) i 


1350 


MOTOR EXCITATION : 


20V, 0-Teak, Squarewave, kOO II z , 2 Phase 


STEADY STATE POWER ( WATTS ) : 


2.5 « 1 350 RT M 


TACHOMETER TYPE i 


P.M. Rotor, 2 Phase Wound Stator, Sinusoidal 
Output, 1 Cycle/Ret/Phase , Amplitude carles 
with speed. 


QUANTITY RUILT ! 


k? 


COMMENTS : 


AC Squirrel Cage Induction Motor 


PLIGHT HISTORY : 


More than 20 of these units have flown. 
Units with operating time greater than 5 
years are presently functioning In-orblt. 



Figure 32. BendLx Aerospace Reaction Wheel Technical Data. 
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Figure 33. Gravity gradient effect with 3 RVV (X axis) lx= 15, ly = IS, Iz=3 
K, = 2 . 0 , 7 ,= 1.0 
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Figure 34. Gravity gradient effect with 3 RVV (Y axis) Ix= 15, Iy= 18, Iz=3 

K y = 2.0 , t ,= 1.0 
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Figure 35. 



Gravity gradient effect with 3 RW (Z axis) Ix= 15, Iy= 18, Iz=3 

a; = 2.o , t,= i.o 
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Figure 36. 



Gravity gradient effect with 3 RW (X axis) Ix = 95, Iy = 98, lz=3 

K , = 2.0 , t , = 1.0 
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Figure 37. 



Gravity gradient effect with 3 RW (Y axis) lx = 95, ly = 98, lz^3 

a; = 2.0 , Xy — l.o 
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Figure 38. 



Gravity gradient effect with 3 RW (Z axis) lx=95, ly = 98, Iz=3 
K,= 1.0, t,= 1.0 
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Figures 39-41 have moments of inertia of I x — 15, I y — 18, and 4 = 3 with 
K x = K y = .2, K 2 = .1, and r x = T y = x 2 = 1.0. The reduction of K x , K y and K 2 reduces the 
reaction wheel torque generated. 

Figure 39 shows 0(position), 4>{ velocity), and reaction wheel torque (T rwi ) 
with respect to time. Figures 40 and 41 show 6, 6, T rwe , and t {/, \p , T rwy , , respectively. 
Stabilization, within the + 1° criteria, occurs for all axes in under 200 seconds. Com- 
parison of the magnitudes of the torque values of Figures 33-35 and Figures 39-41 for 
the same moment of inertia values shows a decrease of roughly a factor of five. 

Increasing the moments of inertia to I x = 95, I y = 98 and 4 = 3, while hold- 
ing K x = K y = .2, K 2 = .1, and r, = z y = r. = 1.0, produces Figures 42-44. This last group 
of reaction wheel figures shows that s;650 seconds are required to dampen the system 
to the + 1° stabilization criteria. The tradeoffs between the different reaction wheel 
combinations are essentially time response against torque/power. The faster the desired 
response, the greater the magnitude of the restoring torque, and ultimately, the more 
power that will be consumed in achieving this response. 
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Figure 39. Gravity gradient effect with 3 RW (X axis) Ix= 15, 
K, = .2 , t„= 1.0 



'= 18, Iz= 3 
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Figure 40. 



Gravity gradient effect >\ itli 3 RW (Y axis) Ix= 15, Iy = 18, Iz=3 

K y = .2 , Ty= 1.0 
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Figure 41. 



Gravity gradient effect with 3 RVV (Z axis) Ix= 15, Iy = 18, Iz=3 

A >.1 , t ,= 1.0 
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Figure 42. Gravity gradient effect with 3 RW (X axis) Ix = 95, Iy = 98, Iz=3 

a; = . 2 , t,= i.o 
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Figure 43. 



Gravity gradient effect with 3 (Y axis) RYV Ix=95, Iy = 98, Iz=3 

K,= . 2 , r,= 1.0 
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Figure 44. Gra\ity gradient effect with 3 RW (Z axis) Ix = 95, ly = 98, lz=3 

K = -1 ,*,= 1.0 
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d. Effects of Magnetic torquers 

For the magnetic torquing simulations, the initial conditions are 0, 9. and 
0 position errors = 5° , and 0. 6. and 0 = 0. The initial restoring torque is provided 
by the gravity gradient effect combined with the rigid body dynamics. For the following 
simulations, one of the two conditions must exist prior to energizing the magnetic coils 
to generate magnetic torque: 

1) position greater than .5° and velocity greater than 0, or 

2) position less than —.5° and velocity less than 0. For the following figures, 
the magnetic moment varies from 9.0427 to .5 Wb-m. At the maximum value of 
9.0427W-m. N = 400 turns, 1 = .5 amp, and A = nr 2 where r is the radius at .12m. Subse- 
quent values are obtained by reducing the current drawn. 

The magnetic field in component form for an equatorial orbit can be ap- 
proximated by eqs(l 10-112). Since the earth's magnetic field is highly dependent on the 
satellite's position. ORION must have the capability to determine its position in order 
to effectively utilize magnetic torquing. For simulation purposes, the satellite is assumed 
to be directly over Greenwich at time t = 0. The magnetic coils are turned on and off at 
the specified strength as needed. Varying the strengths of the magnetic moments gener- 
ated would enhance performance, but is beyond the scope of this thesis. The assumption 
that magnetic torque may be generated in any direction due to vector addition is not 
always valid; magnetic torque may not be generated in the direction of the magnetic 
field. 

Figures 45-47 utilized magnetic moments of m x — m y — m 2 = 9.0427 with 
moments of inertias I x = 15, I y = 18, and l 2 = 3. For Figure 45 on page 60, the gravity 
gradient restoring torque takes approximately S00 seconds before meeting a condition 
necessary for energizing the magnetic coils. Due to the small moments of inertia, the 0 
restoring torque is needed for only a short time. Figures 46 and 47 show the 6 and 0 
axes essentially stabilized within the ± 1° criteria after the the initial magnetic torque 
pulse. For the 6 axis, this pulse occurs at the 1200 seconds and lasts for approximately 
300 seconds. The 9 axis has a lower frequency response and, therefore, takes longer to 
to meet the condition necessary to energize the magnetic coils. Due to its relatively small 
moment of inertia, the 0 axis chatters slightly before stabilizing. A smaller magnitude 
of restoring torque eliminates this problem. 

Increasing the moments of inertia to l x = 95, l y = 98, and /, = 3 while leav- 
ing the magnetic moments m x , m y , and m 2 constant at 9.0427wb-m, produces Figures 
48-50. As expected, more restoring torque is needed to stabilized the 0 and 9 axes forc- 
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ing the magnetic coils to remain energized for a longer period of time. The time needed 
for all axes to be within the ± 1° criteria is limited by the d axis. Stabilization occurs 
within 2S00 seconds while for the case with the smaller moments of inertias, stabilization 
occurs at approximately 2100 seconds. 



59 





rMMocci is iirt tGMitmi 11-15 11*10 t r- 9 . •• 




7 



mag—x 



• » 



> it* 

1 

O >»*l 

M *.« 1.1 1.1 M t.l l.t M 1.1 l.l l.l 

Hr* HI* SCI 

tl IS Met IGMI0OWI 



Figure 45. Gravity gradient with magnetic torquing (X axis): Magnetic 

moment = 9.0427W-m lx = 15, ly = 18, Iz=3. 
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Figure 46. Gravity gradient with magnetic torquing (Y axis): Magnetic 

moment = 9.0427W-m Ix= 15, ly= 18, lz = 3. 
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Figure 48. Gravity gradient with magnetic torquing (X axis): Magnetic 

moment = 9.0427W-m Ix = 95, I y = 98, Iz= 3. 
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Figure 49. Gravity gradient with magnetic torquing (Y axis): Magnetic 

moment = 9. 0427W-m lx = 95, Iy=98, Iz = 3. 
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Figure 50. Gravity gradient with magnetic torquing (Z axis): Magnetic 

moment = 9. 0427W-m lx =95, Iy = 98, Iz= 3. 
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Reducing the strength of the magnetic moments reduces the torque that is 
generated and, thereby, slows the response time to stabilization. For Figures 51-56, the 
magnetic moments m x = m y = m 2 = 1.0. Figures 51-53 have moment of inertias I x = 15, 
4= IS , and 4 = 3. Stabilization for the 6 axis occurs within approximately 3600 sec- 
onds. Increasing the moments of inertias to I x = 95, I y = 98, and I z = 3 produces Figures 
54-56. Stabilization within the + 1° criteria for the 6 axis occurs at approximately 1 1000 
seconds. The length of time that the magnetic torqers are energized increases as the 
strength of the moment decreases. Figures 57-62 show the effects of magnetic torquing 
using a maximum magnetic moment of .5\Vb-m for each coil. Figures 57-59 show that 
for low moments of inertia (I, = 15, I y = IS, and 4 = 3), stabilization occurs within the 
± 1° within 4800 seconds. Increasing the moments of inertia to 4 = 95, 4 = 98, and 
4=3 holding the maximum magnetic moment = .5Wb-m, produces Figures 60-62. In 
Figure 61, the 6 axis may take a prohibitively long time. 
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Figure 51. Gravity gradient with magnetic torquing (X axis) 
moment = l.OW-m lx = 15, ly= 18, Iz=3. 
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Figure 52. Gravity gradient with magnetic torquing (Y axis): Magnetic 

moment = l.OW-m Ix = 15, Iy= 18, Iz=3. 
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Figure 53. Gravity gradient with magnetic torquing (Z axis): Magnetic 

moment = l.OW-m Ix= 15, Iy= 18, Iz= 3. 
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Figure 54. Gravity gradient with magnetic torquing (X axis): Magnetic 

moment = 1.0 W-m lx = 95, 1 y = 98, lz= 3. 
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Figure 55. Gravity gradient with magnetic torquing (Y axis): Magnetic 

moment = l.OW-m Ix = 95, I y = 98, Iz = 3. 
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Figure 58. Gravity gradient with magnetic torquing (Y axis): Magnetic 

moment = .50W-m lx= 15, Iy= 18, Iz=3. 
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Figure 59. Gravity gradient with magnetic torquing (Z axis): Magnetic 

moment = .50W-m Ix= 15, ly= 18, Iz=3. 
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Figure 60. Gravity gradient >\ ith magnetic torquing (X AXIS): Magnetic 

moment = .50 W-m lx = 95, Iy = 98, Iz=3. 
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Figure 61. Gravity gradient with magnetic torquing (Y AXIS): Magnetic 

moment = .50W-m Ix=95, I y = 98, Iz=3. 
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Figure 62. Gravity gradient with magnetic torquing (Z AXIS): Magnetic 

moment = .50W-m Ix=95, Iy = 98, Iz= 3. 
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3. Gravity Gradient capture (from tumbling mode) 

Increasing the inertia will reduce the tumbling rate by conforming to the con- 
servation of angular momentum, I 0 co 0 = /,&>,. As the inertia increases, a corresponding 
decrease in angular rate must be realized. This is useful in obtaining the initial gravity 
gradient capture from a tumbling mode. Figure 63 on page 80 shows that for a slow 
tumbling satellite, the change in the moment of inertia by extending the boom may be 
sufficient to allow gravity gradient capture. The satellite is initially tumbling at a rate 
of 5 rotations per orbit about the pitch axis. The extension of the boom increases the 
moment of inertia about the y axis by a factor of 6; therefore, by the conservation of 
angular momentum, the angular rate must decrease by a factor of 6. This reduces the 
rotations per orbit to approximately 5/6 per orbit well below the one rotation per orbit 
rate necessary to acheive gravity gradient capture. Figure 63a shows the initial tumbling 
rate of 5 rotations per orbit. Figure 63b shows the corresponding decrease in the tum- 
bling rate as the moment of inertia is increased. Figure 63c shows the increase in mo- 
ment of inertia as the boom is extended with respect to time. 

When an initial tumbling rate is higher than the rate that can be captured solely 
through the extension of the boom, magnetic torquing may aid in gravity gradient cap- 
ture. By creating torque in the direction opposite to that of the angular velocity, mag- 
netic torquing can be used to slow the rate of the tumbling vehicle. 
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Figure 63. Gravity gradient capture by extension of boom 
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V. CHAPTER 5 



A. CONCLUSIONS 

1. Gravity gradient stabilization (only) 

Gravity gradient stabilization alone is not a viable attitude control method for 
ORION. This is due to gravity gradient efiect providing little yaw restoring torque. In 
addition, gravity gradient stabilization provides little dampening efiect on the roll and 
pitch axes. When augmented by other control schemes, gravity gradient stabilization 
provides sufficient passive restoring torque to warrant its use. 

2. Gravity gradient with 3 reaction wheels 

a. Advantages 

Gravity gradient stabilization augmented by three reaction wheels, each 
aligned along the body's axes, provides excellent attitude control results both in time 
response and accuracy. 

b. Disadvantages 

Power consumption may be prohibitively high for this method. At steady 
state, each reaction may use up to 6 watts while drawing considerably more during the 
acceleration or deceleration phase. If a steady state drain of 18 watts (6 for each re- 
action wheel) is acceptable, then this method should be considered. The projected end- 
of-life power budget for ORION is 60 watts [Ref. 1, p. 132]. Using this figure as an 
estimate, a steady state draw of 18 watts leaves 42 watts available for mission support. 
From Figure 3 on page 3. an average of 34 watts is required for mission support. The 
size and weight of three reaction wheels may be prohibitively high as well. Another 
disadvantage is that thrusters must be used to remove momentum when the reaction 
wheel(s) become saturated. 

3. Control moment gyro (CMG) 

a. Advantages 

Although not discussed in detail, control moment gyros are capable of 
generating large restoring torques. 

b. Disadvantages 

The weight requirements for a control moment gyro are prohibitively high 
when compared to reaction wheels when relatively small torques are invoved. As a small 
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multi-purpose satellite, ORION does not need the magnitude of restoring torque neces- 
sary to warrant the use of control moment gyros. 

4. Gravity gradient with three magnetic torquers 

a. Advantages 

Three othogonally oriented magnetic torquers have the advantage of gen- 
erating restoring torque in almost any direction based on vector addition. In this man- 
ner, the corrective torque may be produced to counter any disturbance from any 
direction as long as it is small enough not to tumble the vehicle. Another advantage for 
magnetic torquing is that the torquers only draw power while providing restoring torque. 
This amounts to a considerable power savings when compared to the reaction wheel. 
The power drawn by the magnetic torquers during operations is based on the amperes 
used. This allows for flexibility of design since doubling the number of turns and halving 
the current drawn still provides the same magnitude of the magnetic moment. 

b. Disadvantages 



Magnetic torquers are altitude dependent. The magnitude of the earth's 
magnetic field is a function of At high altitudes, magnetic torquing is not effective. 

/V 

In addition, magnetic torquing is not always available in all directions. Recalling that 
magnetic torque is T mat = M x B, magnetic torque may not be generated in the direction 
of the B field due to vector cross product relationships. Due to the continuously chang- 
ing B field with respect to the satellite, this is a minor disadvantage. The desired restoring 
torque can soon be generated as soon as the B field direction changes. 

5. Thrusters 

a. Advantages 



At high altitudes, thrusters and momentum exchange devices are the only 
control methods available for accurate stabilization. In addition, thrusters are the most 
effective method for stabilizing a vehicle with a high tumbling rate. 
b. Disadvantages 

Thrusters are necessary for satellites; however, any fuel used for attitude 
control or momentum removal is less fuel available for orbit transfer and maintainance. 
For a small, lightweight satellite with limited fuel storage capacity, fuel resources must 
be economized. 

6. Overview 

In the final analysis, the attitude control system employed will be based on 
specific mission requirements. Only when the ORION mission is known can proper 
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weights be given to each of the attitude control system design factors, most noticeably, 
power consumption, pointing accuracy, altitude, and desired lifetime. 
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APPENDIX GRADIENT STABILIZATION PROGRAM 



TITLE 3 AXIS STABILIZATION WITH GRAVITY GRADIENT EFFECT 






* DEFINITION OF TERMS 



* ALT=ALTITUDE OF SATELLITE (METERS) 

* R1=RADIUS OF EARTH (METERS) 

* R=RADIUS OF EARTH + SATELLITE ALT (METERS) 

* RO=DENSITY OF AIR AT THE ORBITAL ALTITUDE 

* WO=ANGULAR RATE OF SATELLITE WITH RESPECT TO EARTH 

* D=DRAG FORCE=l. 7E+14*(RO/R) 

* PSI, THETA, PHI ARE ANGLES ABOUT THE X,Y,Z, AXES (RADIANS) 

* PSIC ,THETAC,PHIC ARE ANGLES ABOUT THE X,Y,Z AXES (DEGREES) 

* PSID ,TD , PHID ARE ANGULAR VELOCITIES ABOUT THE X,Y,Z AXES (RADS/SEC) 

* (Bl) ,B( 2) ,B( 3) ARE THE ANGULAR ACCELERATION ABOUT THE 

* X,Y,Z AXES (RADS/SEC**2) 

* IX, IY , IZ ARE THE MOMENTS OF INERTIA ABOUT THE X,Y,Z AXES (KG-M**2) 

* TX , TY , TZ ARE THE SUMMATION OF THE DISTURBANCE TORQUES AND 

* RESTORING TORQUES EXCLUDING GRAVITY GRADIENT EFFECTS 

* TMX , TMY , TMZ ARE MAGNETIC TORQUES ABOUT THE X,Y,Z AXES 

* TRWX , TRWY , TRWZ ARE REACTION WHEEL TORQUES ABOUT THE X,Y,Z AXES 

* TSX ,TSY,TSZ ARE SOLAR PRESSURE TORQUES ABOUT THE X,Y,Z AXES 

* TAX,TAY,TAZ ARE AERODYNAMIC TORQUES ABOUT THE X,Y,Z AXES 

* LX,LY,LZ IS THE DISTANCE BETWEEN THE CENTER OF PRESSURE AND 

* CENTER OF MASS (M) 

* MX ,MY,MZ=SPACECRAFT' S GENERATED MAGNETIC MOMENT WITH RESPECT TO 

* THE X , Y , Z AXES (WB-M) M=NIA WHERE M=MAGNETIC MOMENT, 

* N=# OF TURNS, I=CURRENT (AMPS), AND A=AREA (M**2). FOR 

* THE MAXIMUM VALUE OF M=9.0427, N=400, I=. 5A,RADIUS=. 12M. 

* MG=STRENGTH OF EARTH'S MAGNETIC DIPOLE=8. 0E15 WB-M AT R=0 

* BX,BY,BZ=EARTH'S MAG FIELD WITH RESPECT TO X,Y, Z(WB/M**2) 

* LAMS=RIGHT ASCENSION FROM GREENWICH @ TIME=0, IN DEGREES ( LONGITUDE) 



ARRAY B( 3) 

CONST ALT=1000E+03,RO=7. 0E-10,R1=6378. E+03,. . . 

IX=95. , IY=98. , IZ=03. 00 , LAMS0=0,... 

MX=9. 0478 , MY=9. 0427 , MZ=9.0427, MG=8. OE+15 , . . . 
KX=2. , KY=2. , KZ=1. , TAU1=1. , TAU2=1. , TAU3=1. ,... 
TSX=0. , TSY=0. , TSZ=0. , TAX=0. , TAY=0. , TAZ=0. 



Vf Vf '«V 

* PHI, THETA, AND PSI ARE INIT 5 DEGREES = .08725 RADIANS OF ERROR* 

J- J- j. J- J- y. J- J- J- j- j. j. j . j. j. j. j. j, j. j. ^ j. J, J- j- Jj. J- J- J- J. J- ^ J. J. J. y . J- ^ J. J. J- J. .U Jj. 



INCON ITHETA=. 08725, IPHI=. 08725, IPSI=. 08725, ILAMS=0,PHIC=5. 
THETAC=5 . , PSIC=5. 

* 



DERIVATIVE 

R=R1+ALT 

D=l. 7E+14*(RO/R) 
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WO=SQRT( 3. 987E+14/(R**3) ) 



NOSORT 

LAMS=LAMSO+WO*TIME*( 57 . 3) 
10 IF (LAMS. GE. 360) THEN 
LAMS=LAMS-360. 

GO TO 10 
ENDIF 

* 



* 

****&**************************&****&****************** 

* THIS SECTION ADDS THE REACTION WHEEL EQUATIONS * 

* * * * * * * * * * * * * * * * * * * * * * * * * * * * ‘V * * * * * * * * ‘V * * * * * * * “V * “V * * * * * * * 



* TRWX=-KX*(TAU1*PHID + PHI) 

* TRW Y =-KY*(TAU2 *TD + THETA) 

* TRWZ=-KZ*(TAU3*PSID + PSI) 

***********jrjsjD REACTION WHEE L* *** *** * *** ** *** * ****** * * * 
******************************************************* 



* THIS SECTION ADDS THE MAGNETIC TORQUE EQUATIONS * 

* AND SOLVES AN EXAMPLE FOR THE EARTH'S MAGNETIC * 

* FIELD (EQUATORIAL ORBIT ASSUMED) * 



* * * * * * * * * * * * * * * 



* ** * * * * * 



******************************** 



BX=(MG/(R**3) )*SIN( . 2042)*SIN( ( ( LAMS*PI )/ 180) -5. 0789) 

BY=-(MG/(R**3))*coS(. 2042) 

BZ=-(MG/(R**3))*SIN(. 2042)*COS( ( ( LAMS*PI )/180) -5. 0789) 



TMAGX=( MZ*BY ) - ( MY*B Z ) 
TMAGY=(MX*BZ)-(MZ*BX) 
TMAGZ=( MY*BX) -( MX*BY) 



* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 

* CONDITIONS FOR TORQUE GENERATION MET? 
****************************** 



IF( PHIC. GT. . 5) THEN 
IF(PHID. GT. 0) THEN 
TMX= - ( AB S ( TMAGX ) ) 
ELSE 
TMX=0 
ENDIF 
ENDIF 

IF(PHIC. LT. -.5) THEN 
IF(PHID. LT. 0) THEN 
TMX=ABS( TMAGX) 

ELSE 

TMX=0 

ENDIF 

ENDIF 

IF(THETAC. GT. . 5) THEN 
IF(TD. GT. 0) THEN 
TMY=-(ABS(TMAGY) ) 
ELSE 
TMY=0 
ENDIF 
ENDIF 

IF(THETAC. LT. 5) THEN 
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IF(TD. LT. 0) THEN 
TMY=AB S ( TMAGY ) 

ELSE 
TMY=0 
END IF 
END IF 

IF(PSIC. GT. . 5) THEN 
IF(PSID. GT. 0) THEN 
TMZ=-( ABS(TMAGZ) ) 

ELSE 

TMZ=0 

ENDIF 

ENDIF 

IF( PSIC. LT. -.5) THEN 
IF(PSID. LT. 0) THEN 
TMZ=ABS(TMAGZ) 

ELSE 

TMZ=0 

ENDIF 

ENDIF 

*********************£ ND MAGNETIC TORQING SECTI ON* * * * * * * * * * * ********* 



TX=TMX+TRWX+TSX+TAX 
TY=TM Y+TRWY+TS Y+TA Y 
TZ=TMZ+TRWZ+TSZ+TAZ 



*******FOR GRAVITY GRADIENT EFFECT ONLY SET TX=TY=TZ=0************** 



B(1)=-4*(W0**2)*((IY-IX)/IZ)*PHI+W0*((IX+IZ-IY)/IX)*PSID+(TX/IX) 
B ( 2 ) = - 3*( WO** 2) *( ( IX - I Z ) / I Y) *THETA+TY/ I Y 

B(3)=-(W0**2)*((IY-IX)/IZ)*PSI-W0*((IX+IZ-IY)/IZ)*PHID+(TZ/IZ) 



PHID=INTGRL( 0. ,B(1)) 
TD=INTGRL( 0. ,B( 2) ) 
PSID=INTGRL( 0. ,B(3)) 



CPHI=INTGRL( 0. ,PHID) 
CTHETA=INTGRL( 0. ,TD) 
CPSI=INTGRL( 0. , PSID) 



PHI=IPHI+CPHI 

THETA=ITHETA+CTHETA 

PSI=IPSI+CPSI 



*******' 






********************************* 



*PHI, THETA, PSI ARE = TO INITIAL CONDITIONS AS ANNOTATED WITH AN "i" 
-'PRECEDING VARIABLE+CALCULATED VALUES AS ANNOTATED WITH THE "C" IN FRONT 

* * * * V * * * * * * * * * * * * * * * * * * ‘V * * * * * * * * -5V * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 



* 



PHIC=PHI*57. 3 
THETAC=THETA*5 7 . 3 
PSIC=PSI*57. 3 



****************************************************************** 
* USE FOLLOWING IF GRAVITY GRADIENT EFFECT ONLY * 
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VcVf*yr*^V*******Vf*V«V**VcVf*Vf*Vf**Vc****V«Vyr****Vc'VfVf******yc’Vc*VfVf**Vc*VeVf*V<'**Vc , *Vf 

*ONTRL FINTIM=20000,DELT=200. 

*RINT 200. , PHIC, THETAC, PSIC, PHIDC, TDC , PSIDC 
*AVE 200., PHIC, THETAC, PSIC, PHIDC, TDC, PSIDC 

*RAPH (G1,DE=TEK618) TIHE(UN=SEC) ,PHIC(UN=DEGREES) 

*ABEL (G1,DE=TEK618) PHI(DEG) VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G2 ,DE=TEK618) TIME(UN=SEC) , THETAC (UN=DEGREES) 

*ABEL (G2,DE=TEK618) THETA(DEG) VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G3,DE=TEK618) TIME(UN=SEC) ,PSIC(UN=DEGREES) 

*ABEL (G3,DE=TEK618) PSI(DEG) VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G4,DE=TEK618) TIME(UN=SEC) ,PHIDC(UN=' DEG/SEC 1 ) 

*ABEL ( G4 , DE=TEK6 18 ) PHID VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G5,DE=TEK618) TIME(UN=SEC) ,TDC(UN=' DEG/SEC ' ) 

*ABEL (G5,DE=TEK618) TDC VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH ( G6 ,DE=TEK618 ) TIME(UN=SEC) , PSIDC (UN=' DEG/ SEC’ ) 

*ABEL (G6,DE=TEK618) PSIDC VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*************£^0 GRAVITY GRADIENT ONLY GRAPHS ******************* 

* 

yc , ****yf*yf*******ycyf***ycyf***^v*^v************‘V**yc*yr*yr^v***yf* , 5V'5V* , 5V , 5V 

* USE THE FOLLOWING IF REACTION WHEELS ARE USED * 

* ‘k *>V “V Vr * * * Vr ‘V •>( * * ‘V Vc * * Vc * Vc Vf * Vc ‘V * Vr * ‘Wr V* •k Vc k * * k *>V k k kk k kk * * kk k kk k *Y ?We * 

*ONTRL FINTIM=400 ,DELT=5. 

*RINT 100. , PHIC, THETAC, PSIC, TX, TY, TZ 

*AVE 5. , PHIC, THETAC, PSIC, TX, TY, TZ, PHID, TD, PSID 

*RAPH (G1,DE=TEK618) TIME(UN=SEC) ,PHIC(UN=DEGREES) 

*ABEL (G1,DE=TEK618) PH I ( DEG) VS TIME (GGJilOOOK) IX=95 IY=98 IZ=3. 00 
*RAPH (G2,DE=TEK618) TIME(UN=SEC) ,THETAC(UN=DEGREES) 

*ABEL (G2,DE=TEK618) THETA(DEG) VS TIME (GG^IOOOK) IX=95 IY=98 IZ=3. 00 
*RAPH (G3 ,DE=TEK618) TIME(UN=SEC) ,PSIC(UN=DEGREES) 

*ABEL (G3 ,DE=TEK618) PSI(DEG) VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G4,DE=TEK618) TIME(UN=SEC) ,TX(UN=NM) 

*ABEL (G4,DE=TEK618) TXRW VS TIME (GG?1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G5,DE=TEK618) TIME(UN=SEC) ,TY(UN=NM) 

*ABEL (G5 ,DE=TEK618) TYRW VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G6,DE=TEK618) TIME(UN=SEC) ,TZ(UN=NM) 

"ABEL (G6,DE=TEK618) TZRW VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G7 ,DE=TEK618) TIME(UN=SEC) ,PHID(UN=* RAD/SEC* ) 

*ABEL ( G7 , DE=TEK6 18 ) PHID VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH ( G8 , DE=TEK6 18 ) TIME(UN=SEC) ,TD(UN=’RAD/SEC’ ) 

*ABEL (G8 ,DE=TEK618) TD VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G9,DE=TEK618) TIME ( UN=SEC ), PS ID(UN=' RAD/SEC* ) 

*ABEL (G9 ,DE=TEK618) PSID VS TIME (GG^IOOOK) IX=95 IY=98 IZ=3. 00 

*********************end reaction wheel graphs**^****************** 

/V k 



* USE FOLLOWING GRAPHS IF MAGNETIC TORQUING IS UTILIZED * 

*ONTRL FINTIM=4000. ,DELT=5. 

*RINT 100., PHIC, THETAC, PSIC, WO,B( 1) ,B(2) ,B(3) ,LAMS 

*AVE 5. , PHIC, THETAC, PSIC, LAMS, TX, TY,TZ,BX, BY, BZ, PHID, TD, PSID 

*RAPH (G1,DE=TEK618) TIME(UN=SEC) ,PHIC(UN=DEGREES) 

*ABEL (G1,DE=TEK618) PHI(DEG) VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G2,DE=TEK618) TIME ( UN=SEC ) , THETAC ( UN=DEGREES ) 

*ABEL (G2,DE=TEK618) THETA(DEG) VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G3,DE=TEK618) TIME ( UN=SEC ), PSIC (UN=DEGREES) 

*ABEL (G3,DE=TEK618) PSI(DEG) VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
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*RAPH (G4,DE=TEK618) TIME(UN=SEC) ,LAMS(UN=DEGREES) 

*ABEL (G4,DE=TEK618) LAMS(DEG) VS TIME (GG@1000K) IX=95 IY=98 IZ=3.00 
*RAPH (G5 ,DE=TEK618) TIME(UN=SEC) ,TX(UN=N-M) 

*ABEL (G5,DE=TEK618) TX VS TIME (GG@1000K) 

*RAPH (G6 ,DE=TEK618) TIME(UN=SEC) ,TY(UN=N-M) 

*ABEL (G6,DE=TEK618) TY VS TIME (GG@1000K) 

*RAPH (G7,DE=TEK618) TIME(UN=SEC) ,TZ(UN=N-M) 

*ABEL (G7,DE=TEK618) TZ VS TIME (GG(aiOOOK) 

*RAPH (G8,DE=TEK618) TIME(UN=SEC) ,BX(UN=W-M) 

*ABEL (G8,DE=TEK618) BX VS TIME (GG(aiOOOK) 

*RAPH (G9 ,DE=TEK618) TIME(UN=SEC) ,BY(UN=W-M) 

*ABEL (G9,DE=TEK618) BY VS TIME (GG(aiOOOK) 

*RAPH (G10,DE=TEK618) TIME(UN=SEC) ,BZ(UN=W-M) 

*ABEL (G10 ,DE=TEK618) BZ VS TIME (GG(aiOOOK) 

*RAPH (Gil ,DE=TEK618) TIME(UN=SEC) ,PHID(UN='RAD/SEC' ) 

*ABEL (Gil ,DE=TEK618) PHID(RADSEC) VS TIME (GG(aiOOOK) IX=95 IY=98 IZ=3. 00 
*RAPH ( G 1 2 , DE=TEK6 18) TIME(UN=SEC) ,TD(UN=' RAD/SEC ' ) 

*ABEL (G12 ,DE=TEK618) TD(RADSEC) VS TIME (GG<aiOOOK) IX=95 IY=98 IZ=3. 00 
*RAPH (G13 ,DE=TEK618) TIME(UN=SEC) ,PSID(UN=’ RAD/SEC ' ) 

*ABEL (G13 ,DE=TEK618) PSID(RADSEC) VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 

END 

STOP 
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